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Abstract 

In this paper, we investigate effects of the minimal length on quantum tunnelling from spherically 
symmetric black holes using the Hamilton-Jacobi method incorporating the minimal length. We 
first derive the deformed Hamilton-Jacobi equations for scalars and fermions, both of which have the 
same expressions. The minimal length correction to the Hawking temperature is found to depend 
on the black hole’s mass and the mass and angular momentum of emitted particles. Finally, we 
calculate a Schwarzschild black hole’s luminosity and find the black hole evaporates to zero mass 
in infinite time. 
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INTRODUCTION 


The classical theory of black holes predicts that anything, including light, couldn’t escape 
from the black holes. However, Stephen Hawking first showed that quantum effects could 
allow black holes to emit particles. The formula of Hawking temperature was first given 
in the frame of quantum field theoryfl]. After that, various methods for deriving Hawking 
radiation have been proposed. Among them is a semiclassical method of modeling Hawking 
radiation as a tunneling effect proposed by Kraus and Wilczekj^, 3], which is known as the 
null geodesic method. Later, the tunneling behaviors of particles were investigated using 


the Hamilton-Jacobi method 


4-6]. Using the null geodesic method and the Hamilton-Jacobi 


rT 


method, much fruit has been achieved[]7H18]. The key point of the Hamilton-Jacobi method 
is using WKB approximation to calculate the imaginary part of the action for the tunneling 
process. 


On the other hand, various theories of quantum gravity, such as string theory 


quantum gravity and quantum geometry, predict the existence of a minimal length 


19 


oop 


21 ], 


The generalized uncertainty principle(GUP) 22] is a simply way to realize this minimal 
length. 


An effective model of the GUP in one dimensional quantum mechanics is given by 23), 24] 

P 


Lfk(p ) = tanh 


M 


( 1 ) 
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LfOj(E) = tanh ( - 

\M f 

where the generators of the translations in space and time are the wave vector k and the 
frequency c o, Lf is the minimal length, and LfMf = h. The quantization in position 
representation x = x leads to 

k = —id x , co = id t . (3) 

3 

Therefore, the low energy limit p -C Mf including order of gives 



P = , 

E = m f 1 “ JP, d ‘ ) ' 


(4) 

(5) 


where we neglect the factor |. From eqn. (II]) , it is noted that although one can increase 
p arbitrarily, k has an upper bound which is j-j. The upper bound on k implies that that 
particles could not possess arbitrarily small Compton wavelengths A = 2ir/k and that there 
exists a minimal length ~ Lf. Furthermore, the deformed Klein-Gordon/Dirac equations 
incorporating eqn. (J4]) and eqn. (J5]) have already be obtained in 23], which will be briefly 
reviewed in section mi So [23| provides a way to incorporate the minimal length with special 
relativity, a good starting point for studying Hawking radiation as tunnelling effect. 

The black hole is a suitable venue to discuss the effects of quantum gravity. Incorporating 


GUP into black holes has been discussed in a lot of papers[25 


The thermodynamics of 


black holes has also been investigated in the framework of GUP |29!, [30(]. In 31], a new form 
of GUP is introduced 


o 

o 

(6) 

p l = k l (l — ak + 2ct 2 k 2 ) , 

(7) 


where p a is the modified four momentum, k a is the usual four momentum and a is a small 
parameter. The modified velocity of photons, two-dimensional Klein-Gordon equation the 
emission spectrum due to the Unruh effect is obtained there. Recently, the GUP deformed 
Hamilton-Jacobi equation for fermions in curved spacetime have been introduced and the 
corrected Hawking temperatures have been derived 3214361]. The authors consider the GUP 
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of form 


%0i, 

Pi = p 0i (1 + Pp 2 ), 


( 8 ) 

(9) 


where Xo i, and poi satisfy the canonical commutation relations. Fermions’ tunnelling, black 
hole thermodynamics and the remnants are discussed there. 

In this paper, we investigate scalars and fermions tunneling across the horizons of black 
holes using the deformed Hamilton-Jacobi method which incorporates the minimal length 
via eqn. (J3J) and eqn. (J5]) . Our calculation shows that the quantum gravity correction is 
related not only to the black hole’s mass but also to the mass and angular momentum of 
emitted particles. 

The organization of this paper is as follows. In section ttU from the modified fundamental 
commutation relation, we generalize the Hamilton-Jacobi in curved spacetime. In section 
IIII1 incorporating GUP, we investigate the tunnelling of particles in the black holes. In 
section El we investigate how a Schwarzschild black hole evaporates in our model. Section 
El is devoted to our conclusion. In this paper, we take Geometrized units c — G — 1, where 
the Planck constant h is square of the Planck Mass m p . We also assume that the emitted 
particles are neutral. 

II. DEFORMED HAMILTON-JACOBI EQUATIONS 

To be generic, we will consider a spherically symmetric background metric of the form 



where where / (r) has a simple zero at r = r/,, with f (r/,) being finite and nonzero. The 
vanishing of / (r) at point r = Th indicates the presence of an event horizon. In this section, 


we will first derive the deformed Klcin-Gordon/Dirac equations in flat spacetime and then 
generalize them to the curved spacetime with the metric (HUD . 

In the (3 + 1) dimensional flat spacetime, the relations between (pi,E) and (ki,u) can 
simply be generalized to 



( 11 ) 


( 12 ) 
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where, in the spherical coordinates, one has for k 


, J 9 d ip d 
k = —i | r— + + 


( 13 ) 


dr r 80 r sin 9 d(j) 

Expanding eqn. dill) and eqn. (fT2|) for small arguments to the third order gives the energy 
and momentum operator in position representation 

K 2 


E = m 1 1 - M~f‘ 


_ h 
p = - 

i 


r d r 


h 2 dl 

Mf 


+ »|* 


Vdl \ 

r 3 M 2 ) 


d<h 


h 2 dl 


r sin 8 r 3 sin 3 9M 2 


(14) 

(15) 


where we also omit the factor |. Substituting the above energy and momentum operators 
into the energy-momentum relation, the deformed Klein-Gordon equation satisfied by the 
scalar field with the mass m is 

E 2 (f) = p 2 (j) + m 2 (j), (16) 

where p 2 = p ■ p. Making the ansatz for (j) 

'if 


= exp 


h 


(17) 


and substituting it into eqn. (TT6|) . one expands eqn. (TT6l) in powers of h and then finds that 
the lowest order gives the deformed scalar Hamilton-Jacobi equation in the flat spacetime 


W? 1 + 


2 (d t lf 

M f 


— ( d r I ) 2 ( 1 + 


2 (d r I) 2 \ (delf 


r 2 sin 2 9 


1 + 


Mf 

2 ( d^iy 

r 2 sin 2 8M 2 f 


1 + 


2 ( d e lf 

r 2 M 2 


= m 


(18) 


which is truncated at O 


Mi 


Similarly, the deformed Dirac equation for a spin-1/2 fermion with the mass m takes the 
form as 

( 70 -E + 7 • p - m) = 0, (19) 

where { 70 , 70 } = 2, {7 a ,7&} = —2 5 a b, and { 70 , 7 a} = 0 with the Latin index a running 
over r, 9 ,and 0. Multiplying ( 70 -E + 7 • p + m) by eqn. (TT9|) and using the gamma matrices 
anticommutation relations, the deformed Dirac equation can be written as 


E 2 0 = (p 2 + m 2 )f - ^ ! a idb\ papb'ip . 


( 20 ) 
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To obtain the Hamilton-Jacobi equation for the fermion, the ansatz for ^ takes the form of 


V’ = exp ( ^ ) v, 


( 21 ) 


where v is a vector function of the spacetime. Substituting eqn. (12T|) into eqn. (120]) and 
noting that the second term on RHS of eqn. (1201) does not contribute to the lowest order of 
h, we find the deformed Hamilton-Jacobi equation for a fermion is the same as the deformed 
one for a scalar with the same mass, namely eqn. (ITSlh Note that one can use the deformed 
Maxwell’s equations obtained in 23] to get the deformed Hamilton-Jacobi equation for a 
vector boson. However, for simplicity we just stop here. 

In order to generalize the deformed Hamilton-Jacobi equation, eqn. (fTSlh to the curved 


spacetime with the metric (ITU , we first consider the Hamilton-Jacobi equation without GUP 
modifications. In ref. we show that the unmodified Hamilton-Jacobi equation in curved 
spacetime with ds 2 = g^dx^dx" is 


g^d^ldul — m 2 = 0 . 

Therefore, the unmodified Hamilton-Jacobi equation in the metric (fTOl) becomes 


o djy 


f(r) {d r If 


(■delf Wf 


= m 


( 22 ) 


(23) 


f (r) J K v ' ' r 2 r 2 sin 2 9 

On the other hand, the unmodified Hamilton-Jacobi equation in flat spacetime can be ob¬ 
tained from eqn. (jT8|) by taking Mf —> 00 , 

(<delf (<V ) 2 


0 dtlf ~ (d r lf ~ 


= m 


(24) 


r 2 r 2 sin 2 9 

Comparing eqn. (12 3 p with eqn. (12T|h one finds that the Hamilton-Jacobi equation in the 
metric (TT01) can be obtained from the one in flat spacetime by making replacements d r I —» 
y/f (r)d r I and d t I -^= in the no GUP modifications scenario. Therefore, by making 
replacements d r I —» \Jf (r)d r I and d t I —> the deformed Hamilton-Jacobi equation 

in flat spacetime, eqn. (TT8jh leads to the deformed Hamilton-Jacobi equation in the metric 
(□HP , which is to O 


M, 


{dtlf 

f( r ) 


1 + 


2 (djy 

f (r) M 2 


/ (r) ( d r lf 1 + 


2 f{r) (drlf 
M f 


(delf 


1 + 


2 (d e lf 

r 2 M 2 f 


{d*rf 


r 2 sin 2 9 


1 + 


2 (V) 2 

r 2 sin 2 9M ? 


= m 2 . 


(25) 
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III. QUANTUM TUNNELLING 


In this section, we investigate the particles’ tunneling at the event horizon r = r/, of 
the metric m where GUP is taken into account. Since the metric (TTUT) does not depend 
on t and 0, d t and 8$ are killing vectors. Taking into account the Killing vectors of the 
background spacetime, we can employ the following ansatz for the action 


I = -ut + W (r, 6) + 


(26) 


where c o and p$ are constants and they are the energy and the ^-component of angular 
momentum of emitted particles, respectively. Inserting eqn. (l26j) into eqn. (l25]h we hnd 
that the deformed Hamilton-Jacobi equation becomes 


pUi + 2f(r)p ' 


Mj 


where p r = d r W, pe = 8qW and 


f 2 (r) 


oj 2 \1 + 


2cn 2 


/ (r) M2 


- / (r) (m 2 + A) 


(27) 


2 P 2 e \ Pi 


2 P 2 d 


\ _ Po_ | i I | I j <p [ i I 

r 2 \ r 2 M| J r 2 sin 2 6 y r 2 sin 2 OMj 


Since the magnitude of the angular momentum of the particle L can be expressed in terms 
of p e and p#, 


one can rewrite A as 


Solving eqn. d2TD for p r to O () gives 


d r W T = 


1 2 

“ Pe + • 2 a ’ 
sm 6 

(28) 

L 2 / 1 \ 

(29) 

7? + 0 (mj)- 


'f(r)\ 


LU 2 1T 


2w 2 


/ (r) M 


/ (r) (m 2 + A) x 


1 - 


\ M 2 /(r) _ 


cu 2 1 + 


2u 2 


f (r) M 2 


/ (r) (m 2 + A) 


(30) 


where +/— represent the outgoing/ingoing solutions. In order to get the imaginary part of 
W±, we need to hnd residue of the RHS of eqn. (j30jf at r = rh by expanding the RHS in a 
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Laurent series with respect to r at r = r^. We then rewrite eqn. (1301) as 


d r W T = 




"h/ M + ^2 


/ 2 (r) (m 2 + A) x 


\ 


/2 (r) “ M? 




/ 2 (r) (m 2 + A) 


( 31 ) 


Using / (r) = /' ( 77 J (r — 77 J + - ^ (r — r/j ) 2 + (9 ((r — r^) 3 ), one can single out the ^ 
term of the Laurent series 

±a i 


—rh 


d r W T ~ 


r - 


where we have 


a; 


a_i = 


/'Oh) 


1 + 1 ” + “ 


O 


M* 


(32) 


(33) 


Using the residue theory for semi circles, we obtain for the imaginary part of W± to O 


Im W± = ±— 


7 TCO 


f (r h ) 


1 + M} [m + “ 


(34) 


However, when one tries to calculate the tunneling rate T from Im W±, there is so 


called “factor-two problem ’ 1 


ing a “temporal contribution’' 


38 [. One way to solve the ’’factor two problem” is introduc- 


14 




39 


40]. To consider an invariance under canonical 


transformations, we also follow the recent work 144161. l39l . 140 ] and adopt the expression 
(j)p r dr = f p+dr — f p~dr for the spatial contribution to T. The spatial and temporal 
contributions to T are given as follows. 

Spatial Contribution: The spatial part contribution comes from the imaginary part of 
W (r). Thus, the spatial part contribution is proportional to 


exp 

= exp 


— — Im (j) p r dr 


= exp 


— — Im 

h 

27 tuj 


f (Th) 


P'1 dr — / p r dr 




(35) 


Temporal Contribution: As pointed ill Ref. jld 3d iii the temporal part contribn- 
tion comes from the ’’rotation” which connects the interior region and the exterior region of 


the black hole. Thus, the imaginary contribution from the temporal part when crossing the 

























































horizon is Im ( ujAt out ’ m ) = where /c = - ^ is the surface gravity at the event horizon. 
Then the total temporal contribution for a round trip is 

2ttoj 


Im( " A() = /'(r,)' 

Therefore, the tunnelling rate of the particle crossing the horizon is 


(36) 


T oc exp 


= exp 


f Im (a;At) + Im (j p r dr 


4:71UJ 

f (r h ) 


1 + 


M) 




m 2 + 


(37) 


This is the expression of Boltzmann factor with an effective temperature 

r K) h 


T = 


47T 


(38) 


1 + Bj( m2 + f)’ 

where T 0 = is the original Hawking temperature. 

For the standard Hawing radiation, all particles very close to the horizon are effectively 
massless on account of infinite blueshift. Thus, the conformal invariance of the horizon make 
Hawing temperatures of all particles the same. The mass, angular momentum and identity of 
the particles are only relevant when they escape the potential barrier. However, if quantum 
gravity effects are considered, behaviors of particles near the horizon could be different. For 
example, if we send a wave packet which is governed by a subluminal dispersion relation, 
backwards in time toward the horizon, it reaches a minimum distance of approach, then 
reverse direction and propagate back away from the horizon, instead of getting unlimited 


blueshift toward the horizon 


41 


42], Thus, quantum gravity effects might make fermions 


and scalars experience different (effective) Hawking temperatures. However, our result shows 
that in our model, the tunnelling rates of fermions and scalars depend on their masses and 


angular momentums, but independent of the identities of the particles, to O 


Mi 


In other 


words, effective Hawking temperatures of fermions and scalars are the same to O 
our model as long as their masses, energies and angular momentums are the same. 


Mj 


m 


IV. THERMODYNAMICS OF BLACK HOLES 

For simplicity, we consider the Schwarzschild metric with / (r) = 1 — with the black 
hole’s mass, M. The event horizon of the Schwarzschild black hole is ry = 2 M. In this 
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section, we work with massless particles. Near the horizon of the the black hole, angular 
momentum of the particle L ~ prh ~ c orh- Thus, one can rewrite T 

To 


T 


1 + 


2ui 2 ’ 


(39) 


where Tn = 


8nM 


for the Schwarzschild black hole. As reported in 43], the authors obtained 


the relation uj > between the energy of a particle and its position uncertainty in the 
framework of GUP. Near the horizon of the the Schwarzschild black hole, the position un¬ 


certainty of a particle will be of the order of the Schwarzschild radius of the black hole 44] 
5x ~ rh- Thus, one finds for T 

T -(40) 


1 + 


m p 

2 M 2 M'j 


where we use h = m 2 . Using the first law of the black hole thermodynamics, we find the 
corrected black hole entropy is 

dM 


S = 


! T 

A iTTml 


4m 2 


M f 


In 


A 

167T 


(41) 


where A = 47rr^ = 167 tM 2 is the area of the horizon. The logarithmic term in eqn. (]4T]i 
is the well known correction from quantum gravity to the classical Bekenstein-Hawking 
entropy, which have appeared in different studies of GUP modified thermodynamics of black 


holes 


27 


29 


45 


51]. In general, the entropy for the Schwarzschild black hole of mass M in 


four spacetime dimensions can be written in form of 


S — -f a In 

4 


A 

167T 


+ o 


MJ 

A 

M ? 


2 M 2 f M 2 \ 

where a = in our paper. Neglecting the terms O i—fj 
three scenarios depending on the sign of a. 


m eqn. 


(42) 

, there could be 


1. cr = 0 : This case is just the standard Hawking radiation. The black holes evaporate 
completely in finite time. 

2. a < 0 : The entropy S as function of mass develops a minimum at some value of 
Mmin . This predicts the existence of black hole remnants. Furthermore, the black 
holes stop evaporating in finite time. This is what happened in 
is consistent with the existence of a minimal length. 


3-29 


45 


51] , which 
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3. a > 0 : This is a subtle case. In the remaining of the section, we will investigate how 
the black holes evaporates in our model. 


For particles emitted in a wave mode labelled by energy u and L, we find from eqn. (15711 


that 


52| 


(Probability for a black hole to emit a particle in this mode) 

= exp ^ x (Probability for a black hole to absorb a particle in the same mode), 


where T is given by eqn. (hi8[) . Neglecting back-reaction, detailed balance condition requires 
that the ratio of the probability of having N particles in a particular mode with u and L to 
the probability of having N — 1 particles in the same mode is exp (—. One then follows 
the standard textbook procedure to get the average number n^.L in the mode 


n u ,L = n (^) , 


where we define 


n (x) = 


exp x — (— 


and e = 0 for bosons and e = 1 for fermions. In 


y 


(43) 


(44) 


53], counting the number of modes per 


frequency interval with periodic boundary conditions in a large container around the black 
hole, Page related the expected number emitted per mode n Uj L to the average emission rate 
per frequency interval dn £ L by 

rl/n T rh ,? 

(45) 


dn u x du 

= n *,L- 


dt ’ 2irh 

Following the Page’s argument, we find that in our model 


dn. 




dt 




du; dp r 
dp r 2irh 




du 
2i rh’ 


(46) 


where is the radial velocity of the particle and the number of modes between the wavevec- 
tor interval ( p r ,p r + dp r ) is where p r = d r I is the radial wavevector. Since each particle 
carries off the energy u, the total luminosity is obtained from dn '“ t L by multiplying by the 
energy u and summing up over all energy u and L , 

du 


L — ^2 ( 2 ^ + 

i=o ^ 


un Ut i 


2i rh’ 


(47) 


where L 2 = (l + 1) lh 2 and the degeneracy for l is (21 + 1). However, some of the radiation 
emitted by the horizon might not be able to reach the asymptotic region. We need to 
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consider the greybody factor |T) (a;)| 2 , where T t (a;) represents the transmission coefficient of 
the black hole barrier which in general can depend on the energy oj and angular momentum 
l of the particle. Taking the greybody factor into account, we find for the total luminosity 


L 


^2(21 + 1) ( \Ti (w)| 2 U7l Uyl -^r. 
1=0 


(48) 


Usually, one needs to solve the exact wave equations for \T t (cu)| 2 , which is very complicated. 
On the other hand, one can use the geometric optics approximation to estimate |T) (a;)| 2 . 
In the geometric optics approximation, we assume a; 3> M and high energy waves will be 
absorbed unless they are aimed away from the black hole. Hence we have |7j(c<;)| 2 = 1 
for all the classically allowed energy uj and angular momentum l of the particle. In this 
roximation, the Schwarzschild black hole is just like a black sphere of radius R = 3 3 / 2 M 
, which puts an upper bound on /(/ + !) h 2 , 



l (l + 1) h 2 ^ 27MV. 


(49) 


Note that we neglect possible modifications from GUP to the bound since we are interested 
in the GUP effects near the horizon. Thus, the luminosity is 

c27M 2 uj 2 


L = 


Lodoj 
2nh 3 


T*M 2 r°° 

2nh 3 In 


n 

r*27 


U 


77T U + 


l l(l + i)h 2 
M 2 2 M 2 


d[l(l + 1) h 2 ] 


u 3 du / n [u (l + a 2 ux)] dx, 


(50) 


where we define u=%,y= l -^f and a = ^ = For M » we have 

a <C 1 and hence the luminosity is 


27 


-Tq A / u 3 n (u) du, 


(51) 


327 v 2 h 3 

which is just the Stefan’s law for black holes. Therefore for large black holes, they evaporate 

2 

in almost the same way as in Case 1 until M ~ , when the term a 2 ux starts to 

dominate in eqn. (1501) . Then the luminosity is approximated by 


T 0 4 M 2 f°° 


-27 


2-Kh 3 


u 3 du / n ( a 2 u 2 x ) dx 


2 M 2 M A 


f 


r*27 


nm. 


6 


^dv / n ( v 2 x ) dx, 


(52) 
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where v = au. Not worrying about exact numerical factors, one has for the evaporation rate 


dM 

dt 


= —L ~ —A 


M 2 

M]’ 


(53) 


where A > 0 is a constant. Solving eqn. 


for M gives M 


M 


A. The evaporation 
The black hole then 


rate considerably slows down when black holes’ mass M ~ 8 J^ P M • 
evaporates to zero mass in infinite time, ffowever, the GUP predicts the existence of a 
minimal length. It would make much more sense if there are black hole remnants in the 
GUP models. How can we reconcile the contradiction? When we write down the deformed 


Hamilton-Jacobi equation, eqn. ([25]) . we neglect terms higher than O 


However, when 


M ~ 8 J^ P M f ’ our e ff ec ti ye approach starts breaking down since contributions from higher 


order terms become the same important as these from terms O 


Mi 


. Thus, one then has to 


include these higher order contributions. For example, if there are higher order corrections 
to eqn. (j53j) as in 


dM 

dt 


-A 


M 2 

U 


B 


M 3 

w 


(54) 


where B > 0, one can easily see that there exists a minimum mass M m \ n Mf \ / 4 for black 


l f' 

holes. In another word, the O terms used in our paper are not sufficient enough to 

produce black hole remnants predicted by the GUP. To do so, one needs to resort to higher 
order terms if the full theory is not available. It is also noted that when at late stage of a 
black hole with a > 1, eqn. (139|) becomes 

T 0 


T 


M, 


(55) 


2 M 2 Mj 


which means the tempreture of the black hole goes to zero as the mass goes to zero. 

The GUP is closely related to noncommutative geometry. In fact, when the GUP is 
investigated in more than one dimension, a noncommutative geometric generalization of 
position space always appears naturally 22]. On the other hand, quantum black hole physics 
has been studied in the noncommutative geometry 55]. In 56|, a noncommutative black 
hole’s entropy received a logarithmic correction with a < 0. However, j|37, 58] showed that 
the corrections to a noncommutative schwarzschild black hole’s entropy might not involve 
any logarithmic terms. In either case, the tempreture of the noncommutative schwarzschild 
black hole reaches zero in finite time with remnants left. 


13 














V. CONCLUSION 


In this paper, incorporating effects of the minimal length, we derived the deformed 
Hamilton-Jacobi Equations for both scalars and fermions in curved spacetime based on 
the modified fundamental commutation relations. We investigated the particles’ tunneling 
in the background of a spherically symmetric black holes. In this spacetime configurations, 
we showed that the corrected Hawking temperature is not only determined by the properties 
of the black holes, but also dependent on the angular momentum and mass of the emitted 
particles. Finally, we studied how a Schwarzschild black hole evaporates in our model. We 
found the black hole evaporates to zero mass in infinite time. 
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